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ABSTRACT 

The events recorded by ARGO-YBJ in more than hve years of data collection 
have been analyzed to determine the diffuse gamma-ray emission in the Galactic 
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plane at Galactic longitudes 25° < I < 100° and Galactic latitudes |6| < 5°. The 
energy range covered by this analysis, from ~350 GeV to ~2 TeV, allows the 
connection of the region explored by Fermi with the multi-TeV measurements 
carried out by Milagro. Our analysis has been focused on two selected regions 
of the Galactic plane, i.e., 40° < I < 100° and 65° < I < 85° (the Gygnus 
region), where Milagro observed an excess with respect to the predictions of 
current models. Great care has been taken in order to mask the most intense 
gamma-ray sources, including the TeV counterpart of the Gygnus cocoon recently 
identihed by ARGO-YBJ, and to remove residual contributions. The ARGO-YBJ 
results do not show any excess at sub-TeV energies corresponding to the excess 
found by Milagro, and are consistent with the predictions of the Fermi model 
for the diffuse Galactic emission. From the measured energy distribution we 
derive spectral indices and the differential flux at 1 TeV of the diffuse gamma- 
ray emission in the sky regions investigated. 

Subject headings: cosmic rays-diffuse radiation:Galaxy-disk: methods:observational 


Introduction 


Diffuse gamma-rays are the sum of contributions from several components: the truly 
diffuse Galactic gamma-rays produced by the interaction of cosmic rays, the extragalactic 
background and the contribution from undetected and faint Galactic gamma-ray sources. On 
the Galactic plane the Galactic gamma-rays dominate the other components. The processes 
leading to this emission are the interaction of cosmic nuclei with the interstellar gas through 
the production and decay of secondary 7r° mesons, the bremsstrahlung of high-energy cosmic 
electrons, and their inverse Gompton scattering on low-energy interstellar radiation helds. 
The spectrum of this radiation may provide insight into the propagation and conhnement 
in the Galaxy of the parent cosmic rays, their source distribution and their spectrum at 
the acceleration sites. Since gamma-rays are not deflected by magnetic helds, the diffuse 
component traces cosmic rays and the interstellar environment in distant parts of the Galaxy. 
In addit ion, the Galactic diffuse emission represents the natural background to many different 
signals flMoskalenko et al.ll2004l) . The knowledge of the diffuse emission is necessary for the 
accurate detection of gamma-ray sources, either as point-like or extended. High-quality data 
on the d iffuse emission in the Ga lactic Genter region are needed to constrain the dark matter 
models (IMacias fc Gordonll201411 . Electrons and positrons from astrophysical sources invoked 
to explain the PAMELA results also induce gamma-rays when propagating in the Galaxy. 
The induced gamma-rays represent an additional contribution to the diffuse gamma-ray 
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background fjZhang et alJl2010l ) . 


Galactic diffuse gamma-rays with ene rgies from 100 MeV to a few GeV w ere firstly de¬ 
tecte d by the space-borne detectors SAS-2 flKniffen fc Fichtellll98ll) and GOS-B fjStrong et ah 
19871 ). which revealed the noticeable correlation between the flux of gamma-rays and the 
density of the interstellar medium. Then GOMPTEL and EGRET on board the Comp¬ 
ton Gamma Ray Observatory provided the hrst maps unv e iling the spec trum of the diffuse 
emission from 1 MeV up to 10 GeV (iHnnter et al.l 119971 : IStrongi 120111) . In the GeV en¬ 
ergy range, the EGRET data show in all directions a signihcant excess (the so-called “GeV 
excess”) with respect to the predictions obtained assuming a cosmic-ray flux as that mea¬ 
sured at Earth. This result stimulated intense theoretical studies and many possible ex- 
planations have be en proposed, including harder cosmic-ray spectra throughout the Galaxy 
f Strong et ah 200o[). large contrib ution from high-energy electrons via inverse Gompton scat¬ 
tering (IPorter fc Protheroel 119971 ). or pr oduction of ga mma-rays from annihilation of dark 
matter particles with mass 40—50 GeV IjBi et al.l 120081) . At intermediate Galactic latitudes 
this excess has not been conhrmed by the Larg e Area Telescope ( LAT) on board the Fermi 
Gamma Ray Space Telescope launched in 2008 (lAbdo et al.ll2009al) . With a sensitivity more 
than one order of magnitude better than EGRET, Fermi-LAT mapped the gamma-ray sky 
up to a few hundreds of Ge V with unprecedented accuracy. Observation and discoyery of 


gamma-ray 


discr ete sources flNolan et al.l l2012l: lAckermann et ai.l l2013l: lAcero et al.l 12013 


Abdo et al.ll2013l) as well as the study of the GeV diffuse emission at different Galactic lati¬ 
tudes flAckermann et aLll2012a|) have been the primary targets of these analyzes. At higher 
energies, due to their low fluxes, the diffuse gamma-rays can be efficiently studied only 
by ground-based detectors with large effective areas. The Imaging Atmospheric Gherenkov 
Telesc opes (lAGTs), such as Whipple and HEGRA , have set upper limits in very narrow re¬ 
gions (ILeBohec et al.ll2000l: lAharonian et al.ll20011) . The H.E.S.S. telescope array has carried 
out a survey of the Galactic plane covering the region —75° < / < 60° in longitude at an 
energy thres hold of about 250 GeV, a nd has presented the latitude prohle of diffuse emission 
at |6| < 2° flAbramowski et al.l 120141) . This represents the hrst observational assessment of 
diffuse TeV gamma-rays by a Gherenkov telescope. The spatial correlation of gamma-rays 
with giant molecular clouds in the inner Galaxy has been also investigated, the hardness 
of gamma-ray spectrum in dicating that these gam ma-rays originate from protons and nu¬ 
clei rather than electrons (lAharonian et al.ll2006al) . However, though lAGTs are the most 
sensitive detectors operating in the held of gamma-ray astronomy, limited by their held of 
view (FOV) they are not well suited to observe dihuse gamma-rays that have a large-scale 
structure. On the other hand, air shower arrays providing a large FOV and a very high duty 
cycle look more adequate for sky survey purposes. Many air shower exp eriments have set 


upper limits to dihuse gamma-rays, such as Tibet ASy at TeV energies flAmenomori et ah 
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20061) , and EAS-TOP (Aglietta et al.lll992[) . KASCADE ( Haungs et al.ll2006l) . and CASA- 


MIA flBorione et al.l 119981) at energies >100 TeV. The Milagro detector has made the hrst 
positive observation of the diffuse gamma-ray flux from the Galac tic plane, measurin g the 
integral flux above 3.5 TeV in the region 40° < I < 100°, \b\ < 5° fjAtkins et ah 2005 ). As¬ 
suming that the contribution from inverse Compton scattering is negligible at TeV energies, 
the Milagro measurement is many times higher than expected, indicating the existence of 
a “TeV e xcess” possibly connect ed to the “GeV excess” of diffuse gamma-rays observed by 
EGRET (jProdanovic et al.ll2007l ). Most notably, in a following analysis of seven-years data, 
the Milagro collaboration reported a clear excess of diffuse gamma-rays at 15 TeV median 
energy from the galactic region 65° < I < 85°, s uggesting the presence of active cosmic- 
ray sources accelerating hadrons flAbdo et al.ll2008l) . This interpretation is supported by the 
fact that this longitude interval harbors the Cygnus X star-forming region, rich with possible 
cosmic-ray acceleration sites. However, the Fermi-LAT data do not conhrm at GeV energies 
the broadly distributed ex cess of diffuse emission observed by Milagro at multi-TeV energies 
(lAckermann et ahll2012bl) . Instead, they have found a brig ht extended gamma-ray source 


centered at / = 79°.6 ± 0°.3, the so-called “Gygnus cocoon” flAckermann et ah 


TeV c ounterpart has been recently identihed by the ARGO-YBJ experiment (iBartoli et ah 


201 ih. whose 


2ni4al ). These results suggest that the gamma-ray excess of the cocoon is likely due to a 


population of freshly accelerated hadronic cosmic rays, even though a leptonic or mixed 
origin cannot be discarded. High-energy protons lastly diffusing from this source and inter¬ 
acting with the local gas could generate gamma-rays that contribute to the diffuse excess 
observed by Milagro. To better clarify the interpretation of these results, we have used the 
data collected by ARGO-YBJ in 5.3 yr to measure the diffuse gamma-ray emission in the 
25° < I < 100°, \b\ < 5° Galactic region, then selecting the two regions 40° < I < 100° and 
65° < / < 85° in order to compare the results with the Milagro measurements. Thanks to its 
high-altitude location and its particular layout, ARGO-YBJ is able to operate at an energy 
threshold of about 300 GeV, in such a way providing results that bridge the Fermi GeV 
energies and the multi-TeV range explored by Milagro. 


The plan of the paper is as follows. In Section 2 we describe the ARGO-YBJ detector 
and its performance. In Section 3 we present data selection and data analysis methods. The 
results of the analysis are presented in Section 4. Section 5 is devoted to a summary of the 
results and to the conclusions. 
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2. The ARGO-YBJ experiment 


The ARGO-YBJ detector, hosted in a building at the YangBaJing Cosmic Ray Obser¬ 
vatory (Tibet, China, 90°31’50” E, 30°06’38” N), 4300 m above sea level, has been designed 
for very high-energy (VHE) gamma-ray astronomy and cosmic-ray observations. It is made 
up of a single layer of resistive plate chambers (RPCs) operated in streamer mode, 2.850 
m X 1.225 m each, organized in a modular conhguration to cover a surface of about 5600 
m^ with an active area of about 93%. The RPCs detect the charged particles in air show¬ 
ers with an efficiency > 98%. To improve the shower reconstruction, other chambers are 
deployed around the central carpet for a total instrumented area of 100 m x 110 m. A 
highly segmented readout is performed by means of 55.6 cm x 61.8 cm external electrodes, 
called ’’pads,” whose fast signals are used for triggering and timing purposes. These pads 
provide the digital readout of the detector up to 22 particles/m^, allowing the count of the 
air shower charged particle s without any signih cant saturation up to primary cosmic-ray 
energies of about 200 TeV flBartoli et al.N2012al) . In order to extend the dynamical range 
to PeV energies each RPC is also equipped with two large size pads (139 cm x 123 cm) 
allowing the colle ction of the total charge developed by the particles hitting the detector 
(lAielli et al.ll20121) . The digital output of each pad is splitted in two signals sent to the logic 
chain that builds the trigger and to the 18,360 multi-hit time-to-digital converters, which are 
routinely calibrated with 0.4 ns accuracy b y means of an off-line method using cosmic-ray 
showers flHe et ahl 120071: lAielli et al.ll2009al). Mo re de tails about the dete ctor and the RPC 


performance can be found in lAielli et ahl (120061) and lAielli et ahl fj2009bl) . The detector is 


connected to two independent acquisition systems correspond ing to two differen t operation 
modes, referred to as the shower mode and the scaler mode (lAielli et al.ll2008l) . The data 
used in this paper were recorded by the digital readout in shower mode. This mode is im¬ 
plemented by means of an inclusive trigger based on the time correlation between the pad 
signals, depending on their relative distance. In this way the data acquisition is triggered 
when at least 20 pads in the central carpet are hred in a time window of 420 ns. By means 
of this trigger the energy threshold for gamma-induced showers can go down to 3 00 G eV 
with an effective area depending on the zenith angle (see Figure 1 in iBartoli et ahl (120131) 1. 


The event reconstruction follows a standard procedure allowing a detailed space-time 
reconstruction of the shower front, including the position of the shower core and the incident 
direction of the primary partic le. A detailed account of the reconstruction algorithm can be 
found in iBartoli et ahl (12011 alibi. l2013l. |2015|). Briefly, the shower core position is obtained 
by htting the lateral density distribution of the shower charged particles with a Nishimura- 
Kamata-Greisen-like function using the maximum likelihood method. The arrival direction 
of the showers is reconstructed by the least squ ares method assuni ing a conical shape for the 
shower front, as described by Equation (1) in (lAielli et al.ll2009al) . which gives the relation 
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between the particle arrival time and the distance to the shower core. The angnlar resolution 
depends on the number of bred pads Npad- The opening angle 'ipio containing 71.5% of the 
events from a point source is about 2° for events with Npad > 20, 1.36° for Npad > 60 and 
0.99° for events with Npad > 100. The number of hit pads Npad is the observable related 
to the primary energy. However, the number of particles at ground level is not a very 
accurate estimator of the primary energy of the single event, due to the large fluctuations 
in the shower development in the atmosphere and to its partial sampling with the limited 
detector area. The primary energy distribution corresponding to different Npad intervals is 
very broad, spanning over more than one order of magnitude for small Npad values. The 
rel ation between N^^and the primary gamma-ray energy of the selected showers is reported 


m 


Bartoli et al.l (120151) . Since the variable Npad does not allow the accurate measurement 


of the primary energy of the single event, the energy spectrum is evaluated by studying the 
global distribution of Npad- The observed distribution is compared with a set of simulated 
ones obtained with different test spectra, in order to hnd out that which better reproduces 
the data (as carried out in Section 4.1). The angular resolution, pointing accuracy, absolute 
energy calibration and detector stability are tested by measuring the cosmic-ray shadow 
cast by the Moon, d etected with a signihcance of 10 standard deviations (s.d.) per month 
( Bartoli et ah 2011b). 


3. Data analysis 


The ARGO-YBJ experiment began taking data in its full conhguration in 2007 Novem¬ 
ber at a trigger rate of 3.5 kHz with a dead time of 4%. It has been operated stably for 
more than hve years, up to 2013 January, with an average duty cycle of 86%, for a total 
effective time of 1670.45 days. For the present analysis, events with zenith angles less than 
50°, corresponding to the declination interval —20° < J < 80°, and Npad > 20 are used. 
A set of standard cuts, applied to the shower core reconstructed position and to the time 
spread of the shower front, have been used to select high-quality data. With this data se¬ 
lection a total of 6.407 x 10^° shower events are observed from the Galactic plane in the 
latitude belt |6| < 15°. The fraction of survived events is about 80%. With these selections 
more back ground cosmic rays than gamma-rays are rejected, implying an increase of the 
sensitivity (iBartoli et al.l 120131) . These data have been used to measure the diffuse emission 
from the regions of the Galactic plane of longitude 25° < / < 100° and 130° < I < 200°. 
Indeed the region 100° < / < 130° is excluded since in the high declination region 6 > 60° 
the Galactic plane runs parallel to the right ascension axis and the contribution from the 
signal could affect the background estimation. The excess of gamma-induced showers is 
obtained following the procedure of the background estimation applied to the ARGO-YBJ 











data as reported in iBartoli et al.l fl2011ai l2013l) . All data are divided into three pad groups, 
20 < Npad < 59, 60 < Npad < 99 and Npad > 100. For each group of pad multiplicity both 
the inner and the outer regions of the Galactic plane are divided into a grid of 0.1° x 0.1° 
bins and hlled with the detected events according to their reconstructed arrival directions 
(event map). The number of cosmic-ray b a.ckground events (bac kground map) is estimated 
by using the direct integration method of iFlevsher et al.l (120041) . The effect of cosmic-ray 
anisotropy on the back ground evaluation ha s been estimated and corrected by applying the 
normalization given in lBartoli et al.l (j2011al ). This procedure is applied to each map bin us¬ 
ing a surrounding region of 16° x 16° in which the estimated background is renormalized to 
the detected events. The ±5° region around the Galactic plane and the 4° x 4°/cos(6) region 
around the Grab Nebula position are excluded from the normalization procedure. However, 
since the diffuse gamma-ray emission extends to more than |&| = 5°, its contribution causes 
an overestimation of the correction related to the cosmic-ray anisotropy. This effect has 
been evaluated using the latitude prohle provided by the Fermi-LAT model for the diffuse 
Galactic emission (see section 4) smeared out with the ARGO-YBJ point spread function 
(PSF). A variation of 15% of this contribution implies on average a variation of the excess 
of about 4%. To investigate the systematic errors related to the extension of the region used 
to evaluate the anisotropy effect, the box size has been varied from 12° x 12° to 20° x 20° 
obtaining an excess variation of ~10%. Both the event and background maps have been 
smoothed with the PSF corresponding to each Npad interval. Then the background map has 
been subtracted to the event map obtaining the event excess map. This map contains events 
from true diffuse gamma-rays as well as from point and extended sources, whereas, due to the 
background subtraction, the isotropic extragalactic emission is canceled out. Indeed, many 
TeV gamma-ray sources of different extension lying on or close to the Galactic plane have 
been detected in the longitude range 25° < I < 100°. Source locations as given in the TeV- 
Cat0are excluded from the analysis. Faint sources (SNR G54.1-I-0.3, VER J2016-I-372, HESS 
J1923-I-141 (W51) and HESS J1943-I-213) have not been masked. Their total contribution 
to the diffuse flux at 1 TeV is estimated to be about 2.5%. 


Taking into account the angular resolution of the detector and the extension of these 
sources, the contribution from a region 4° x 4°/cos(6) centered around each source location 
has been removed. Some boxes include two sources, with the fainter one near its edge 
(as for instance HESS J1849-000 and VER J2019-I-368). Sources distant less than 1°.2, 
as for instance HESS J1857-I-026 and HESS J1858-I-020, have been masked with a unique 
box centered at the median point. The massive star-forming regio n of Gygnus X hosts th e 
extended cocoon, hrstly observed by Fermi-LAT above 1 GeV (lAckermann et al.ll201ll) . 


^http://tevcat.uchicago.edu 
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whose emission at TeV energies has been recently assessed by ARGO-YBJ flBartoli et ah 
2014a|). Due to its extension of about 2°, this region has been masked with a box 6° x 


6°/cos{b) centered on the source position found by ARGO-YBJ. The chosen dimensions of 
these boxes is a compromise between a desired large excluded region, in order to minimize 
the contamination from the sources, and the requirement of not reducing the statistics. With 
this choice the solid angle of the region 25° < I < 100°, |5| < 5° is reduced of about 22%. 
The spillover from these sources outside the masked regions has been estimated by tracking 
their path inside the FOV of ARGO-YBJ. The contamination is calculated bin by bin and 
subtracted from the total excess in the 0°.l x 0°.l bin event excess map. For ARGO J1839- 
0627/HESS J1841-055, ARGO J1907+0627/MGRO J1908+06, and ARGO J20 31+4157 fthe 


Cygnus coco on) the fluxes measured by ARGO-YBJ have been considered flBartoli et ah 


20131. l2014aj). Since the PSF broadens with decreasing energy, this contamination is found 


higher for the first energy bin (corresponding to the group with 20 < Npad < 59), with an 
average value of 14%, while it is 21% in the Cygnus region 65° < I < 85°. 


Results 


The Galactic longitude prohle of the diffuse gamma-ray emission at 600 GeV in the 
latitude belt |6| < 5° is shown in Fig. [Uby the hlled circles. The prohle obtained without 
masking the sources is also plotted as the open circles. Each point represents the hux 
obtained with the spectral analysis reported in Section 4.1 and averaged over 8° longitude 
bins. The hux is evaluated taking into account all the events with Npad > 20. For each group 
of pad multiplicity the event excess measured in each longitude bin is converted to a hux 
using the eh ective areas estimated by means of a full Monte Carlo simu lation of extensive 
air showers (IHeck et al.l Il998h and of the RPC array flGuo et al.l 1201011 . Then a spectral 
analysis is carried out as described in Section 4.1. The negative huxes which appear in the 
outer Galaxy prohle correspond to negative excess values in the event excess map. In this 
case a spectral index -2.7 has been assumed. The signihcance of the excess measured in the 
25° < I < 100° region is 6.9 s.d., while no excess is detected in the outer Galaxy region 
130° < / < 200°. Fig. |2]shows for |5| < 15° the Galactic latitude prohle of the excess in bins 
of 2°. The hlled circles show the results after masking the sources, while the open circles 
show the results without masking. The continuous line in these plots represents the hux 
at 600 GeV provided by the standard Fermi-LAT model for the dihuse Galactic emission 
galJ2yearp7v6ju0.fits (hereafter Fermi-DGE) to which we refer for comparison with the 
ARGO-YBJ data. This model, available at the Fermi Science Support Centei@, has been 


^http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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used to generate four Fermz-LAT So urce Catalogs flNolan et al.ll2012l: lAckermaun et alJl2ni3 


Acero et al.ll2013l: lAbdo et al.ll2013[) . For the first time a diffuse flux measured by a ground- 
based detector overlaps and can be compared with results from direct measurements. In 
the region 25° < I < 100° we hud a satisfactory general agreement between the ARGO-YBJ 
data and the fluxes predicted by Fermi-DGE, mostly in the inner range 40° < I < 90°. The 
maximum deviations <2.5 s.d. are observed at three values outside this interval. In addition 
to statistical fluctuations, systematic uncertainties related to the background evaluation, 
imperfect modeling of the Galactic diffuse emission and other effects as, for instance, an 
energy-dependent diffuse flux from unresolved sources, can contribute to this discrepancy. 
The amount of such uncertainties is addressed in the next section, where the results of a 
spectral analysis of the flux detected in this region are presented. In the following sections, 
the spectral analyses concerning two selected subregions, 40° < I < 100° and 65° < I < 85°, 
and the upper limit to the diffuse flux in the outer Galaxy are reported and discussed. 


4.1. The Galaxy Region 25° < I < 100°, |6| < 5° 


The total number of shower events recorded in this region is 7.92 x 10®. To carry out a 
spectral analysis with a distribution of the nuni b er of events in excess as a function of Npad 
we follow the method described in lAielli et al.l fl2010l) . Sampling events are generated by 
simulations in the energy range from 10 GeV to 100 TeV assuming the spectral index of a 
power law as a parameter and taking into account the detailed ARGO-YBJ detector response. 
The £t to the data is made by comparing the measured excess in each pad multiplicity interval 
with simulations. A differential spectral index -2.80±0.26 is found. The corresponding 
median energies of the events recorded in the Npad intervals 20 < Npad < 59, 60 < Npad < 
99, Npad Y 100 are 390 GeV, 750 GeV, and 1.64 TeV, respectively. Since the median 
energies depend on the spectral index, these values are affected by an uncertainty of about 
30%. This result is shown in Fig. |3] (dots). Upper limits from HEGRA, Whipple, and 
Tibet ASy experiments are also shown. The solid line represents the expectation according 
to the Fermi-DGE model. This model is dehned between 50 MeV and 600 GeV. Above 
10 GeV the spectrum follows a power law with spectral index about -2.6 and has been 
extended with the same slope to TeV energies as a guide for the eye (dashed line). The 
Fermi-DGE model is based on template fits to the all sky gamma-ray data and includes 
an Inverse Gqmpton component generated by the GALPROFo cosmic-ray propagation code 
(jStrong et al.ll2000t IVladimirov et al.l 120111) . It arises from an accurate comparison of data 
to the sum of many contributions including the ones from detected sources. It is not easy 


^ ht tp: / /sourceforge. net /projects/galprop 
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to asse ss the uncertainty associated to the predictions of the model. In lAckermann et ah 


fl2012a() a grid of models is considered, which represents well the gamma-ray sky with an 
agreement within 15% of the data, although various residuals at a ^30 % level are found, 
both at s mall and large scale s . Sim ilar systematic errors are quoted in iMacias fc Gordon 
(I 2 OI 4 J) and iGordon fc MaciasI (120131 ) in an analysis concerning a very narrow region around 
the Galactic center. Above 40 GeV the accuracy of the modeling is limited primarily by the 
photon statistics and the diffuse emission has been derived by extrapolating the emissivities 
measured at lower energies. The average systematic uncertainty at high energies is expected 
to be greater than 10%, which is the systematics affecting the determination of the effective 
area flAckermann et ahll2012a) . Two main systematic uncertainties can affect the ARGO- 
YBJ flux estimate: one on the background and the other on the absolute scale energy. The 
systematics on the background evaluation has been discussed in Section 3. The energy scale 
reliability has been checked by studying the westward shift of the cosmic-ray shadow cast by 
the Moon due to t he geomagnetic held . At TeV energies the total absolute energy scale error 
is less than 13% (iBartoli et al.ll2011bh and the corresponding systematic error on the flux 
normalization would be about 23%. Minor contributions to the systematic error come from 
the uncertainty on the residual contamination of the masked sources, the ARGO-YBJ and 


H.E.S.S. fluxes being q uoted with about 30% systematic uncertainty (IBartoli et al.ll2ni2b 


^u 

Aharonian et al.ll2006bl) . from the uncertainty (<4%) on the detector efficiency and from the 


systematic error of about 5% affecting the effective area estimate. We combine these various 
errors in quadrature to obtain a total systematic error of ~27%. The estimated ARGO-YBJ 


flux at 1 TeV is (6.0 ±1.3) x 10 TeV 


T-l 


— 2 —1 
cm s 


sr 13% lower than the prediction 


based on the FermADGE extrapolation. Taking into account the whole uncertainties, we 
deem the ARGO-YBJ data set consistent with the model predictions. 


4 . 1 . 1 . Treatment of the unresolved sources 


However, it is worth noting that part of the detected signal could originate in faint 
sources that are unresolved because of their low flux. Indeed, the Galactic region that we 
are studying hosts many potential gamma-ray sources, mainly supernova remnants (SNRs) 
and pulsar wind nebulae (PWNe). Shell-type SNRs and SNRs interac ting with molecular 
cloud s form an established source class in VHE gamma-ray astronomy (iHinton V Hofmann 
20091 ). although whether the na ture of their emi s sion is predominantly hadronic or lep- 
tonic is still a matter of debate ( Yuan et ah 2012 : Volk fc Berezhko 2013). The TeV flux 
depends on the energy available for shock acceleration, on the distance and on many en- 
vironmental parameters . The majority of the identihed Galactic TeV sources are PWNe 
(Ide Oha-Wilhelmi et al.l 120131) . Ninety Galactic PWNe and PWN candidates are reported 
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bv iKargaltsev et al.l (120131) . of which 51 with VHE associations or possible VHE counter¬ 
parts. The formation of a pulsar wind is still poorly understood and it is not known which 
pulsars are able to drive PWNe and produce VHE radiation. In a PWN, the source of the 
energy of the injected electrons is the pulsar spin-down luminosity. TeV PWNe detected 
with the current instruments, at a sensitivity of about 2% of the Crab flux, are mainly as¬ 
sociated with young and energetic pulsars with spin-d own power E > 10^^ erg s~^, showing 


a TeV luminosity not signihcantly correlated with E (IKargaltsev et al.ll2013l: iKlepser et ah 


20131 ). On the other hand, as mentioned above, several energetic pulsars with prominent 


X-ray PWNe are not detected at TeV energies, suggesting that environmental factors such 
as, for instance, the local energy density of the ambient photon held or the intensity of the 
magnetic held, are relevant for the evolution of TeV PWNe. The H.E.S.S. catalog, which 
reports the results of the Galactic plane survey in the longitude range —110° < I < 65°, 
lists 71 TeV sources, of which about 35% hrmly associated with PWNe and 21% with SNRs. 
Besides a few massive stellar clusters and some binary systems, a large fraction (about 31%) 
remains unide ntihed, with ambiguou s associations or without any plausible counterpart in 
X-ray or radio (jCarrigan et ah 1120131) . 15 TeV sources and candidate sources are recorded in 
the 25° < I < 65°, \b\ < 3.5° Galactic region common to H.E.S.S. and ARGO-YBJ. They are 
associated with PWNe or are unidentihed, but many of them could be old PWNe, still bright 


Karg£ 

iltsev et al. 

2013: 

Acero et al. 

2013) 

(2005' 

°) reports 321 radio-loud pulsars in t 


The ATNF catalog (vl.50, Manchester et al. 


selection ehects, the number of predicted pulsars in the same Galactic region is 4169, while 
inside 25° < I < 65°, |6| < 5° is 4633@. To obtain this predictio n, we used the Galacto- 
centric distribution (model G’) suggested bv iLorimer et al.l (120061) . Therefore, the number 
of TeV PWNe detected at the threshold of 2% of the Grab flux is a very tiny fraction of 
the predicted number of their radio counterparts. An estimate of the contribution to the 
measured diffuse emission from sources with lower fluxes would require a model of a full syn¬ 
thetic population of P WNe based on H.E . S.S. a n d Fermi data and on the current theories 
of PWN evolution 


V-g., 


Bucciantini et al.l (120111) : iMaver et al.l (120121) : iTorres et ahl (l2014l) ). 


This study is beyond the scope of the present paper. However, in order t o have an indication 
of the effect of undetected sources, we can use the results obtained by iGasanova fc Dingus 
(I 2 OO 8 I) . Taking into account only the sources detected by H.E.S.S. above 6% of the Grab 
flux, these authors determine the number-intensity relation dN/dS oc S~‘^, where S is the 


^http://www.inpi-hd.mpg.de/hfm/HESS/pages/home/sources/ 

® http://www.atnf.csiro.au/research/pulsar/psrcat 

®with 1.4 GHz luminosities above 0.1 mJy kpc^ and beaming towards us 
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source integral flux above 200 GeV, as expected for a uniform density distribution of sources 
in a two-dimensional disk. Extrap olating this relation do wn to 2 mCrab flux and assuming 
an average spectral index of -2.3 flKargaltsev et al.l 120131) . we can estimate the cumulative 
differential flux at 1 TeV as a function of the numbe r N of unresolved sourc es with S in 
the range 0.2%-2% of the Crab flux above 200 GeV f Aharonian et ah 2006b h As an ex¬ 
ample, for V = 50 we obtain a contribution of about 9% to the diffuse flux measured by 
ARGO-YBJ. This contribution depends linearly on N, but may be lower if the source count 
distribution flattens at low fluxes. A similar exercise can be a pplied to SNRs. A study of 
the contribution of unresolved shell-type SNRs is carried out by I Volk fc Berezhkol (120131) to 
explain the apparent excess of Fermi data at GeV energies in the inner Galaxy. We note 
that the number of SNRs predicted in the 25° < I < 65 °, \b\ < 5° is 69, using the r adial 
distribution of the surface density of sh ell SNRs g i ven in ICase Sz Bhattacharval (119981) and 
the Galactic height distribution given in IXu et al.l (120051) . 


However, these considerations do not settle definitively this matter. Indeed, the H.E.S.S. 
survey extends only to I = 65° and in a latitude belt |5| < 3.5°, and the sensitivity of the 
instrument is reduced for extended sources. Twenty-two SNRs and 749 PSRs are predicted 
in the re gion 65° < I < 100°, | 5| < 5°. Apart from the C ygnus region observed firstly by 


HEGRA (lAhmonian et al.ll2002h an d then by the Whipple (iKonopelko et a 


(lAlbert et al.ll2008l) and VERITAS (IWeinstein II2009I: lAliu et ah 


l ongitude interval h as been surveyed only by Milagro (lAtkins et al. 


2013 


20071) . MAGIC 


20141) telescopes, this 


200 


Jand ARGO-YBJ 


(iBartoli et al.ll2013l) . the latter providing a study with a sensitivity of 24% Crab units. The 
H.E.S.S. survey shows clearly a strong decrease of VHE gamma-ray sources moving toward 
the outer Galaxy, however the presence of some isolated sources with fluxes below the ARGO- 
YBJ sensitivity cannot be excluded. A single Crab-like source with 10% of the Crab flux gives 
contribution at 1 TeV. The Milagro collaboration has found multi-TeV emission from 


the direction of two gamma-ray pu lsars (PSR J1928-I-1746 and PSR J2030-I-3641) detected 
by Fermi-LAT (lAbdo et al.l I20141) . corresponding to a total contribution of (7 ± 2)% to 
the diffuse flux at 1 TeV. However, for one of these sources only a much lower upper limit 
has been obtained by VERITAS, and the flux measured by Milagro might include some 
additional diffuse emission. 


The Fermi-LAT catalog of sources above 10 GeV IFHL (lAckermann et al.l 120131) pro¬ 
vides a list of TeV candidates, stating that many of them should be detectable with the 
current generation of ground-based instruments, but no candidate is found lying in the 
Galactic region I = 85° to / = 100°. In this paper the Fermi Collaboration presents a study 
of the source populations above 10 GeV to infer the contribution of the resolved and unre- 
so lved sou r ces to both high-latitude and low-latitude diffuse backgrounds. The method used 


m 


Stroud (120071) . based on EGRET data with predictions for the Fermi-LAT, is adopted. 
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Therefore, the source count distribution for Galactic sources is modeled with a power law 
luminosity function ~ L~^'^ with gi ven limits, and with a source distribution in Galactocen- 
tric distance based on the model of iLorimer et al.l (120061) for the pulsar distribution, taken 
as representative of Galactic sources. An exponential scale height of 500 pc is assumed. The 
simulated differential source count dN/dS is then compared with the observed flux distri¬ 
bution of the IFHL sources. Both the source density and the luminosity range are varied 
to obtain alternative models. Thus they estimate that the contribution of sources below the 
Fermi-LAT detection threshold of 5 x 10“^° photons cm~^ to the observed gamma-ray 
intensity above 10 GeV at low latitudes (|6| < 10°, all longitudes), is about 5%. This result 
cannot be easily scaled to the TeV range in the Galactic region considered here. Never¬ 
theless, it could suggest a small contribution from unresolved sources. Indeed, for source 
emission extending to TeV energies with a spectral index of 2.3, the Fermi threshold flux at 
10 GeV corresponds to a flux at 200 GeV about 1.5 times the threshold used in our study . 
Moreover pulsars, which are not expected to contribute to the TeV flux (lAbdo et al.l 120131) . 
account for about half of the Galactic sources used in the Fermi estimate. A similar result 
has been recently found at 1 GeV considering the Fermi-LAT third source catalog 3FGL 
(lAcero et al.N2015[l . In conclusion, we can assume that while the main contribution from 
discrete sources has been removed, a residual contribution from unresolved sources could 
still affect the measured fluxes. 


The ARGO-YBJ data have been used to study the interval of this region, 40° < I < 100° 
(section 4.2), which is not rich in point or extended sources apart from the Gygnus cocoon, 
and, in a separate analysis, the innermost part, 65° < I < 85° (section 4.3), which includes 
the Gygnus region. 


4.2. The Galaxy Region 40° < I < 100°, |6| < 5° 

The total number of shower events collected in this region is 7.39 x 10®. After masking 
the discrete sources and subtracting the residual contribution, an excess with a statistical 
signihcance of 6.1 s.d. above the background is found. The result of the spectral analysis 
provides the flux at three median energies (350 GeV, 680 GeV, and 1.47 TeV, with uncer¬ 
tainties of about 30%) as shown in Fig. |4](dots). The fluxes measured by ARGO-YBJ below 
1 TeV are ~20% larger than what expected by the Fermi-DGE model, but are consistent 
within the experimental uncertainties. Fitting the whole ARGO-YBJ data with a power law 
we obtain a spectral index of -2.90±0.31 with a predicted flux at 1 TeV of (5.2 ± 1.5) x 10“^® 
TeV“^ cm“^ s“^ sr“^, compatible with the extrapolation of the Fermi-DGE model. 

We have used these data to face the “TeV excess” anomaly associated to the Milagro 
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result concerning this Galactic region. In fact, in this region the Milagro detec tor made the 


first measurement of the diffuse TeV gamma-ray flux from the Galactic plane flAtkins et ah 
20051) . The measured flux above 3.5 TeV is (6.8±1.5±2.2) x 10“^^ cm“^ s“^ sr“^, which, once 


connected to the EGRET data with a power law with differential spectral index -2.6, reveals 
a “TeV excess” in the diffuse gamma-ray spe ctrum, the corresponding flux being 5-10 times 
higher than expected flAharonian et al.l 120081) . In order to explain the enhanced gamma-ray 
flux seen by Milagro, other contribu tions to the true diffuse flux have been envisaged and 
discussed in IProdanovic et al.l fj2007|) . and include a harder cosmic-ray spectrum, additional 
flux from unresolved sources, excess gamma-rays by inverse Gompton scattering and photons 
from dark matter annihilation. Using the Milagro data, we converted this integral flux 
to the differential flux plotted in Fig. 0] (triangle). This flux is only 34% greater than 
the value expected from the Fermi-DGE extrapolation, therefore within the experimental 
uncertainties. Moreover, the Milagro result does not take into account the contributions from 
the Gygnus cocoon and from the overlapping point or extended sources TeV J2032-I-4130, 
VER J2019-I-407, and VER J2016-I-372. Minor contributions co me also from the H.E .S.S. 
sources. Following the ARGO-YBJ analysis of the Gygnus cocoon (iBartoli et al.ll2014al) and 
taking into account the width of the latitude band, we can evaluate the fraction of the total 
flux generated by these sources. We find that the discrepancy between the Milagro results and 
the Fermi-EGE predictions is almost canceled out. The full set of measurements obtained 
with ground-based experiments is in agreement with direct observations by Fermi-EPET 


Accor d ing to these results and taking into account that the Fermi-LPET data flAbdo et al. 


2009al. 120101) do not support the high-intensity diffuse em ission observed by EGRET (the 
EGRET “GeV excess”), likely due to instrumental effects fjStecker et al.l 120081) . we rule out 
the evidence of any “TeV excess” requiring additional sources or particle production processes 
other than those responsible for the production of Galactic cosmic rays. 


4.3. The Gygnus Region 

The statistical significance of the excess found in the Galactic region 65° < I < 85°, 
\b\ < 5° is 6.7 s.d. above the background. After masking the discrete sources and the Gygnus 
cocoon and subtracting the residual contributions, an excess of 4.1 s.d. is left. This direction 
points into our spiral arm at the Gygnus star-forming region hosting a giant molecular cloud 
complex. Located at a distance of about 1.4 kpc, this region is rich in potential cosmic-ray 
accelerators such as Wolf-Rayet stars, OB associations, and SNRs. Given its peculiarity, this 
region has been the target of numerous multi wavelength observations, including the high- 
energy measurements of Fermi-LAT (GeV) and Milagro (TeV). Fermi-LAT data have been 
used to study the region of galactic coordinates 72° < I < 88°, |6| < 15°, where the bright 
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and extended cocoon has been observed. The spectral energy distribution of gamma-ray 
emission is shown in Fig. [5] (hlled stars) and includes different contributions from the diffuse 
emission, point sources, and extended objects. A global model taking into account all the 
components reproduces satisfactorily the experimental data and implies that the cosmic- 
ray flux averaged over the scale of the whole Cygnus region is similar to that of the local 
interstellar space. The expected energy spectrum according to the Fermi-DGE model in the 
same spatial region is shown for comparison (dot-dashed line). 


This region also received considerable attention by Cherenkov telescopes and air shower 
arrays that have discovered VHE emission from point and extended sources. Recently, the 


ARGO -YBJ experiment observed the TeV counterpart of the Fermi cocoon flBartoli et ah 


2014a)). In this paper a short summary of the previous observations at TeV energies is also 


reported. 


Exploiting its wide FOV, t he Milagro telescop e measured the diffusion emission at Galac¬ 
tic longitudes 65° < I < 85° (lAbdo et al.l l2007al) . The first paper repo rts the flux in th e 
latitude band |6| < 3° at a median energy of 12 TeV. In a following paper fjAbdo et al.ll2008l) . 
more data have been added and a more refined analysis is applied evaluating the flux at a 
median energy of 15 TeV in the region with 65° < I < 85°, |6| < 2°, as reported (filled 
triangle) in Fig. [5l The measured flux is twice the predictions based on the GALPROP code 
optimized to reproduce the EGRET data. This excess has been attributed to the interaction 
with the interstellar medium of hard-spectrum cosmic rays generated by local sources. For 
comparison, we show (dashed line) the expected energy spectrum for this region according to 
the Fermi-DGE model, not available at that time. The discrepancy is reduced, the Milagro 
flux being about 75% higher than the Fermi template, but enough to suggest the presence 
of an excess. The flux measured by ARGO-YBJ at median energies 440 GeV, 780 GeV, and 
1.73 TeV (with uncertainties of about 40%) and averaged over the latitude band |&| < 5° 
is shown with dots. The three points can be htted with a power law with spectral index 
-2.65±0.44. The estimated flux at 1 TeV is (6.2 ±1.8) x 10“^° TeV“^ cm“^ s“^ sr“^, resulting 
about 10% lower than the Fermi-EGE extrapolation. These data do not show any excess at 
energies around 1 TeV, corresponding to the excess found by Milagro at an average energy 
of 15 TeV. One possible explanation of this discrepancy is that the contribution of all the 
discrete gamma-ray sources was not completely removed from the Milagro data. Indeed, the 
exclusion of discrete sources is of crucial importance. According to the ARGO-YBJ data, 
the flux at 1 TeV injected by the cocoon is of the same order as the diffuse emission flux. 


An alternative explanation could be considered if the spectrum measured in the Gygnus 
region is compared with that measured in the complementary part of the 25° < I < 100° 
region. Adding the data from the regions 25° < I < 65° and 85° < I < 100°, we found an 
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excess of 5.6 s.d. above the background. The measured spectrum has an index -2.89±0.33, 
while the estimated flux at 1 TeV is (6.0 ± 1.7) x 10“^° TeV“^ cm“^ s“^ sr“^. Thus there 
is an indication that the spectrum of the diffuse emission in the Cygnus region could be 
harder than that in the complementary part of the 25° < I < 100° longitude interval. 
Assuming that there is actually a difference, a plausible explanation is that the region of 
about 500 pc around the Cygnus cocoon is more abundant of cosmic rays accelerated by 
a nearby source, which produces also the TeV emission from the cocoon, whose spectrum 
has not yet been steepened by diffusion (lAharonian fc AtovanI Il996l: iGabici et ahl 120091) . 
These runaway cosmi c rays may diffuse to a characteristic length of a few hundred parsecs 
fjCasanova et al.ll2010l ) and interact with the local gas producing gamma-rays with the same 
spectral shape via 7r° decay. In fact, since the hadronic interactions at multi-TeV energies 
are ba sically scale-invariant an d TeV photons are not attenuated by the interstellar radiation 
helds flMoskalenko et ahl 120061) . the gamma-ray spectrum is expected to mimic the cosmic- 
ray spectrum. This interpretati on assumes a hadron ic origin of the gamma-ray emission 
from the cocoon as discussed in iBartoli et ahl (j2014al) . In this scenario, the region around 
the cocoon is expected to contain a mixture of ordinary background cosmic rays and young 
cosmic rays with a harder spectrum, released hrst from the source, which diffuse fastly and 
reach a distance depending upon many factors as the injection history of the source, the 
diffusion coefficient and the interstellar gas density. Thus, the diffuse gamma-ray emission 
may consist of two distinct components produced by these two cosmic-ray populations which 
have different spatial extension and different spectral shape. The superposition of these two 
components may produce concave spectra at TeV energies flGabici et ahl 120091) . accounting 
for the Milagro result in case there is a residual excess after re moving all the source c ontri- 
butions. An interesting application of these concepts is given in iGasanova et ahl (l2010l) . The 
ARGO-YBJ data do not allow the study of this phenomenology. Accurate measurements 
of the diffuse gamma-ray emission at TeV energies on an angular scale of a few degrees are 
necessary. Future experi ments with high er angular resolu tion and large FOV, such as HAWG 
flAbevsekara et al.l 120131) and LHAASO flGao et al.ll2010l ). are expected to probe the spatial 
distribution of the photon flux at TeV energies, providing a detailed r nap of the diffuse 
gamma-ray emission in this region. They will also beneht of the GTA flActis et al.N201l[) 
survey, expected to detect very faint discrete sources, thus providing information useful to 
separate the genuine diffuse emission. 

The hadronic origin of the gamma-ray emission from the cocoon can be also probed 
searching for an excess o f GeV-TeV neutrinos over the atmospheric neutr ino background 
from the Gygnus cocoon f Tchernin et ah 20131: Gonzalez-Garcia et ah 2014 ). Interestingly, 
the energy spectrum of the light component (protons plus Heliu m nuclei) of the pri mary 
cosmic rays from a few TeV to 700 TeV measured by ARGO-YBJ flBartoli et al.ll2012al) and 
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by the hybrid experiment ARGO-WFCTA flBartoli et al.ll2014b|) follows the same spectral 
shape as that found in the Cygnus region. A precise comparison of the spectrum of young 
cosmic rays, as those supposed in the Cygnus region, with the spectrum of old cosmic rays 
resident in other places of the Galactic plane, could help to determine the distribution of the 
sources of the cosmic rays observed at Earth. 


4.4. Outer Galaxy 

No excess has been measured in the outer Galaxy region 130° < I < 200°, \b\ < 5° 
(after masking the Crab Nebula). Assuming a spectral index -2.7 the median energy of all 
the events with Npad > 20 is 700 GeV. The corresponding upper limit at 99% conhdence 
level (C.L.) results 5.7 x 10“^° TeV“^ cm“^ s“^ sr“^ and is shown in Fig. [6l where the 
limits obtained at higher energies by the Tibet ASy (3 and 10 TeV) and Milagro (15 TeV) 
experiments are also reported. The Fermi-DGE flux and its extrapolation are shown for 
comparison. The ARGO-YBJ upper limit is compatible with the Fermi model, providing 
an useful constraint to the Galactic diffuse emission around 1 TeV. 


5. Summary and conclusions 

We analyzed the data recorded by ARGO-YBJ over more than hve years for a total 
live time of 1670.45 days, with the aim of measuring the diffuse gamma-ray emission at TeV 
energies in the Galactic region visible from the Northern Hemisphere. After the application 
of appropriate selection criteria, 6.407 x 10^° high-quality events are found in the Galactic 
latitude belt |6| < 15°. These events have been used to measure the gamma-ray diffuse 
emission in the two Galactic regions 25° < I < 100°, |6| < 5° and 130° < I < 200°, 
\b\ < 5° accessible to the experiment. Fluxes and spectral indexes measured by ARGO-YBJ 
in these Galactic regions are reported in Table [H The standard ARGO-YBJ procedure for 
background subtraction has been applied, including a suitable approach to correct for cosmic- 
ray anisotropy. Great care has been taken in removing the emission from known gamma-ray 
sources by masking out the brightest of them and subtracting the residual contributions. An 
excess of 6.9 s.d. above the background is observed in the innermost region 25° < I < 100°, 
\b\ < 5°, which has been the target of a detailed analysis since the pioneering Milagro 
observations at multi-TeV energies have shown signihcant deviations from the predictions 
based on conventional models of diffuse gamma-ray emission. As a reference for our results 
we used the recent Fermi model for diffuse emission extrapolating it to the TeV region. 
Firstly, we have studied the region 40° < I < 100°, where a “TeV excess” in the diffuse 
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gamma-ray spectrum has been suggested. The ARGO-YBJ data have been analyzed to 
derive the differential ffux at three median energies around 1 TeV. Fitting these points with 
a power law, we found a spectrum steeper than the Fermi-DGE extrapolation, with index 
-2.90±0.31, however consistent to within 1 s.d.. The large error on the spectral index is 
due to the short lever arm of these data and to the poor statistics affecting the highest 
energy point. The average ffux is compatible with the Fermi-DGE extrapolation within the 
statistical and systematic uncertainties. After subtracting the contribution of the gamma- 
ray sources detected later on and thus not taken into account, also the ffux measured by 
Milagro at 3.5 TeV is compatible with the Fermi-DGE extrapolation. Therefore, we cannot 
confirm the existence of any excess at TeV energies. 


Table 1: Diffuse Gamma-ray emission from the Galactic plane for |6| < 5°. 


1 Intervals 

Significance 

Spectral Index 

Energy (GeV) 

Flux“ 

25° < / < 100° 

6.9 s.d. 

-2.80 ±0.26 

390 

8.06 ± 1.49 

750 

1.64 ± 0.43 

1640 

0.13 ±0.05 

1000^ 

0.60 ±0.13 

40° < / < 100° 

6.1 s.d. 

-2.90 ±0.31 

350 

10.94 ±2.23 

680 

2.00 ±0.60 

1470 

0.14 ±0.08 

lOOO'* 

0.52 ±0.15 

65° < / < 85° 

4.1 s.d. 

-2.65 ±0.44 

440 

5.38 ± 1.70 

780 

1.13 ±0.60 

1730 

0.15 ±0.07 

1000^ 

0.62 ±0.18 

25° < / < 65° & 

85° < / < 100° 

5.6 s.d. 

-2.89 ±0.33 

380 

9.57 ±2.18 

730 

1.96 ±0.59 

1600 

0.12 ±0.07 

1000^ 

0.60 ±0.17 

130° < 1 < 200° 

-0.5 s.d. 

- 

- 

< 5.7“ 


Note. The median energies and the corresponding differential fluxes are reported. The errors 
are only statistical. 

“In units of 10“® TeV“^ cm“^ s~^ sr“^. 

^This entry gives the result of the fit to the three data points. 

“99% G.L. at 700 GeV. 

A specific study of the Gygnus region (65° < / < 85°) is motivated by the Milagro 
results showing, at energies >10 TeV, a strong enhancement of the diffuse ffux with respect 
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to the model predictions, suggesting the existence of powerful young accelerators as sources 
of hard-spectrum cosmic rays. Indeed, following the Fermi-LAT discovery of the Cygnus 
cocoon at GeV energies, a TeV counterpart has been reported by ARGO-YBJ, providing 
an intense flux of TeV photons. The ARGO-YBJ results on diffuse emission around 1 TeV 
do not exhibit any excess when compared to the Fermi data at lower energies, suggesting 
that the tail of the cocoon flux above 10 TeV and other contributions from discrete sources 
not completely removed from data could explain the excess found by Milagro. The ARGO- 
YBJ measurements cover the energy range from about 400 GeV - 2 TeV and follow a power 
law with spectral in dex -2.654:0.44, a va lue very close to that found for TeV emission from 
the Gygnus cocoon flBartoli et al.ll2014a|) . Thus the spectrum appears flatter than the one 
found in the whole region 25° < I < 100° once the Gygnus region is excluded. Indeed, in 
the combined region 25° < I < 65° plus 85° < I < 100° the spectral analysis provides an 
index -2.89±0.33. These measurements are affected by large errors and their difference has 
a marginal statistical signihcance. 

The diffuse gamma-ray flux measured by ARGO-YBJ can provide useful hints to con¬ 
strain models of Ga lactic origin of the high-energy neutrino excess reported by the Ice- 
Gube Gollaboration (lAartsen et al.l 120131 ). The origin of this excess above the atmospheric 
neutrino background is unknown. Scenarios invoking an extragalactic origin are favored 
flAhlers fc HalzenI 120141) . Alternative mo dels envisaging differen t possibilities for the Galac¬ 
tic neutrino sources have been proposed (lAhlers fc Mnrasel 120141) . One of these assumes that 
the TeV-PeV diffuse gamma-ray emission from the Galactic plane and at least part of the 
neutrino IceGube excess are produced via the same mechanism, that is, the interaction of 
cosmic rays with the interstellar medium. With this conjecture a clear connection between 
the gamma-ray and neutrino fluxes from the Galactic plane can be established. The TeV- 
PeV gamma-ray emission should provide a good estimate for the 100 TeV neutrino signal 


along the whole Galactic plane. Recent paper s have addressed this scenario flNeronov et ah 


20141: lAhlers fc Mnrasel 120141 IFox et al.l 120131) . The fluxes measured by ARGO-YBJ (see 
Table 1) add a hrm estimate to the available current data. However, only upper limits are 
known for the high-energy diffuse gamma-ray emission (see Introduction), limiti ng the ca¬ 


pabil ity to draw hrm c onclusions. Future gam ma-ray observatories as LHAASO fiGao et ah 


20101) and HISGORE f Tluczvkont et al.l 120121) will operate with high sensitivity up to PeV 
energies. The whole set of TeV-PeV data and an increased statistics of the neutrino signal 
from a deeper IceGube exposure are expected to conhrm or rule out the model of a neutrino 
diffuse emission from the Galactic plane. 


In conclusion, the ARGO-YBJ results concerning the diffuse emission at TeV energies 
in the 25° < I < 100°, \b\ < 5° Galactic region are in agreement with the extrapolation of the 
Fermz-DGE model, implying that the questions raised by the Milagro observations can be 
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answered by taking into account the emission of TeV photons from the Cygnus cocoon and, 
to a minor extent, from discrete sources. A spectral analysis of the data has been carried 
out, showing an energy spectrum softer than that of the Fermi-DGE model, but consistent 
within 1 s.d.. On the other hand, the TeV flux averaged over the Cygnus region 65° < I < 85° 
shows a marginal evidence of a harder spectrum, indicating the possible presence of young 
cosmic rays coming from a nearby source. Only an upper limit has been set to the diffuse 
emission in the outer Galaxy region 130° < I < 200°, |6| < 5°, but compatible with the 
extrapolation of the Fermi-DGE model. 
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Fig. 1.— Galactic longitude profile of the diffuse gamma-ray emission in the Galactic latitude 
interval |6| < 5° at an energy of 600 GeV as obtained from the ARGO-YBJ data. The hlled 
circles show the results after masking the sources, while the open circles show the results 
without the masking. The solid line represents the value quoted by the Fermi-DGE model 
at the same energy, smeared out with the ARGO-YBJ PSF. 
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Fig. 2.— Galactic latitude profile of the diffuse gamma-ray emission in the Galactic longitude 
interval 25° < I < 100° at an energy of 600 GeV as obtained from the ARGO-YBJ data. 
The hlled circles show the results after masking the sources, while the open circles show the 
results without the masking. The solid line represents the value quoted by the Fermi-DGE 
model at the same energy, smeared out with the ARGO-YBJ PSF. 
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Energy(TeV) 

Fig. 3.— The energy spectrum of the diffuse gamma-ray emission measured by ARGO-YBJ 
in the Galactic region 25° < / < 100°, |5| < 5° (dots). The solid line shows the flux in 
the same region according to the Fermi-DGE model. The short-dashed line represents its 
extension following a power law with spectral index -2.6. The EGRET results (squares) in 
the same Galactic region 25° < I < 100°, |&| < 5° and the upper limits quoted by HEGRA 
(99% G.L., 38° < / < 43°, |&| < 2°), Whipple (99.9% G.L., 38.5° < I < 41.5°, \b\ < 2°) and 
Tibet ASy (99%G.L., 20° < I < 55°, \b\ < 2°) are also shown. 
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Fig. 4.— Energy spectrum of the diffuse gamma-ray emission measured by ARGO-YBJ in 
the Galactic region 40° < / < 100°, |6| < 5° (dots). The solid line shows the flux in the same 
region according to the Fermi-DGE model. The short-dashed line represents its extension 
following a power law with spectral index -2.6. The EGRET results (squares) in the same 
Galactic region 40° < I < 100°, |6| < 5° and the flux measurement by Milagro (triangle) in 
the same region are also shown. 
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Fig. 5.— Energy spectrum of the diffuse gamma-ray emission measured by ARGO-YBJ in 
the Galactic region 65° < I < 85°, |&| < 5° (dots). The solid line shows the flux according 
to the Fermi-DGE model. The short-dashed line represents its extension following a power 
law with spectral index -2.6. The EGRET results (squares) in the same region are also 
shown. The Milagro result (triangle) for the Galactic region 65° < I < 85°, \b\ < 2° is also 
given. The long-dashed line and its extension (short-dashed line) represent the flux in this 
region according to the Fermi-DGE model. The spectral energy distribution of gamma-ray 
emission measured by Fermi-LAT in the Galactic region 72° < I < 88°, |6| < 15° is also 
reported (stars). The flux in the same region expected from the Fermi-DGE model is shown 
as a dot-dashed line. 
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Fig. 6.— The 99% C.L. upper limit at a median energy of 700 GeV as obtained by ARGO- 
YBJ for the Galactic region 130° < / < 200°, |6| < 5°. The solid line shows the flux in 
the same region according to the Fermi-DGE model. The short-dashed line represents its 
extension following a power law with spectral index -2.6. The EGRET results (squares) in 
the same region are also shown. For comparison, the upper limits from the Milagro (95% 
G.L., 136° < I < 216°, l&l < 2°) and Tibet ASy (99% G.L., 140° < I < 225°, \b\ < 2°) 
experiments are also reported. 



